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Noise reduction panel arrangement and method of calibrating such a panel arrangement 
. •>-■ — D 



The present invention relates to a noise reduction arrangement comprising: 
a plurality of actuators for generating secondary noise to reduce primary noise 
5 generated by at least one primary source; 

a plurality of sensors for sensing the total amount of noise resulting from the primary 
noise after being reduced by the secondary noise and for generating a plurality of 
sensor signals; 

control means for controlling the actuators based on the sensor signals, 

10 the distance between the plurality of actuators and the plurality of sensors being first and 
second surfaces is selected to have an optimised reduction in power RP of the total amount of 
noise relative to the primary noise within a predetermined frequency band. 

Such a noise reduction arrangement is known from J. Guo, e.a., "Actively created 
quiet zones by multiple control sources in free space", J. Acoust. Soc. Am. 101 (3), March 

15 1997, pp. 1492-1501. This document discloses an arrangement with a series of secondary 
sources on a first line and a series of error sensors on a second line, the first and second lines 
being parallel. The primary concern of this document is to create large areas of quiet zones. 
The document observes that such a requirement can be satisfied if the error sensors are not in 
the near field of the secondary sources. According to the document, the distance between the 

20 second line with the error sensors and the first line with the secondary sources should be 
greater than or comparable to the mutual distances between the secondary sources. Guo e.a. 
only present a model for this two line arrangement. Moreover, in their model, all secondary 
sources are controlled by the output signals of all error sensors. Implementing such a control 
arrangement results in a complex controller with many connections and which turns out to be 

25 rather slow in many applications. 

S.J. Elliott et al., Interaction Between Multiple Feedforward Active Control Systems, 
IEEE Transactions on Speech and Audio Processing, Vol. 2, No. 4, 1994, pp. 521-530 [1] 
describe a noise reduction system having a panel of actuators arranged in a first plane and a 
plurality of error sensors in a second plane. The first and second planes are parallel to one 

30 another. Elliott et al. present a mathematical model of a decentralised adaptive feedforward 
control system. They also present results of some physical examples in which there are two 
actuators and two error sensors. In these examples, Elliott et al. introduce the mutual distances 
between the error sensors and the actuators as important parameters to derive conditions as to 
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when such a system is stable. In the physical examples given, the distance between the two 
planes is about 0.3 times the distance between the two actuators. Elliott et al. do not disclose 
the presence of an optimum distance between the two planes as a function of the mutual 
distance between actuators. 

X. Qui, e.a., A Comparison of Near-field Acoustic Error Sensing Strategies for the 
Active Control of Harmonic Free Field Sound Radiation, Journal of Sound and Vibration, 
1998, 215(1), pp. 81-103 [2], disclose the results of a study to find the best location of an error 
sensor relative to a primary noise source. However, this study is limited to a harmonic sound 
field radiated by a monopole primary source and by a dipole-like pair of primary sources. In 
both cases the actuator is a monopole radiating at the same frequency as the primary source. 
No plurality of actuators and plurality of error sensors arranged in respective planes are 
disclosed. 

An active high transmission loss panel is disclosed in WOA-94/05005. However, in 
this patent document the actuators and sensors are all located in the same plane. 

The present invention is directed to a noise reduction arrangement having a plurality 
of actuators in a first surface and a plurality of error sensors in a second surface in which the 
reduction of noise is optimised as a function of the distance between the surfaces and in which 
the control means are simplified. The surfaces may be planes, like in the arrangement of 
Elliott et al. [1], but they may also deviate from planes. They may, e.g., be slightly curved. 

Thus, the noise reduction arrangement as defined above is characterised in that 

• the plurality of actuators are located in a first two dimensional array in a first surface; 

• the plurality of sensors are located in a second two dimensional array in a second surface 
arranged substantially parallel to the first surface; 

• the plurality of actuators are sub-divided into a plurality of sub-sets of actuators; 

• the control means comprise a plurality of controllers, each controller being arranged to 
receive sensor signals of a sub-set of said plurality of sensors and arranged to control one 
single sub-set of actuators; and 

• said first and second surfaces are arrange at such a mutual distance that the reduction of 
power RP is within the following range: 

0.9xRP max <RP<RP max 
in which RP max is maximum obtainable reduction in power of the total amount of noise 
relative to the primary noise at an optimum distance between said first and second surfaces as 
established by testing, where both RP and RP max arc expressed in decibel, 
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wherein the plurality of actuators are arranged in rows and columns, mutual distances 
between adjacent columns and mutal distances between adjacent rows being equal to a 
predetermined actuator distance d x , the plurality of sensors being arranged in the same way 
as the plurality of actuators, the distance d between the first and the second surfaces 
5 meeting the following condition: 

0.5xd x <d<d x 
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relati v e to t he primary nois o , where both RP and RP mm ar c expressed in decibel. 

The present invention is based on the insight that a maximum reduction shows up 
in the curve representing the reduction of the total amount of sound power relative to the 
primary noise as a function of the distance between the surfaces and that it is not 
5 necessary to have each actuator controlled by the output signals of each of the sensors. 
The actual optimum distance where the maximum occurs depends on several parameters, 
like the number of actuators, the number of sensors, the ratio between these two numbers, 
the actual arrangement of the actuators and the actual arrangement of the sensors. The 
optimum distance can be established by testing while increasing the distance between the 
10 surfaces from 0, while adjusting a predetermined control parameter (6) to maintain 
stability. 

Preferably, each controller is arranged to receive sensor signals of only those 
sensors which are within a predetermined range from said controller. 

In one of the arrangements, the number of sensors equals the number of actuators 
p 15 and equals the number of controllers, each controller receiving one of the plurality of 
sensor signals as input signal and controlling one of the plurality of the actuators. When, 
p in such an arrangement, the plurality of actuators are arranged in rows and columns 

g mutual distances between adjacent columns and mutual distances between adjacent rows 

H are equal to a predetermined actuator distance d x and the plurality of sensors are arranged 

M 20 in the same way as the plurality of actuators, the distance d between the first and the 
second surfaces preferably meets the following condition: 

0.5 x d x £ d £ d x . 

In one embodiment, the arrangement includes a supervising controller for moni- 
toring long-term behaviour of the arrangement and for modifying control parameters of 
25 the controllers in order to ensure overall stability of the arrangement. 

Hereinafter, the invention will be explained with reference to some drawings. The 
drawings and explanation are only given by way of example and are not intended to limit 
the scope of the present invention. 

Figure la shows a front view of a plate provided with 48 actuators and 221 sensors 
30 in front of the plate; 

Figure lb shows a schematic cross section view of the arrangement according to 
figure la along line IB-IB in figure la; 

Figure lc shows a schematic electronic black box circuitry for controlling the 
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actuators based on the sensor signals generated by the sensors; 

Figure 2 shows sound power curves radiated from a plate without control, with 
global control and local control, respectively; 

Figure 3 shows condition numbers for the curves shown in figure 2; 
5 Figure 4 shows sound power curves as a function of frequency for an arrange- 

ment with 48 actuators and 48 sensors, the distance d between the actuator plane and the 
sensor plane being a parameter; 

Figure 5 shows curves of broadband reduction in sound power for the arrange- 
ment of figure 4 taking into account all frequencies f < c/2d x , with c the speed of sound 
1 0 in air and d x the distance between adjacent actuators; 

Figure 6 shows sound power curves as a function of frequency for an arrange- 
ment with 48 actuators and 221 sensors, the distance d between the actuator plane and 
the sensor plane being a parameter; 

Figure 7 shows curves of broadband reduction in sound power for the arrange- 
1 5 ment of figure 6, taking into account all frequencies f < c/2d x ; 

Figure 8 shows sound power curves as a function of frequency for a global con- 
trol arrangement with 48 actuators and 221 sensors, the distance d between the actuator 
plane and the sensor plane being a parameter; 

Figure 9 shows broad band reduction of sound power according to figure 8, 
20 taking into account all frequencies f < c/2d x ; 

Figure 10 shows sound power curves as a function of frequency for an arrange- 
ment in which the sound produced is reflected by a further plate parallel to the plate sup- 
porting the actuators, the reflection coefficient R being a parameter; 

Figure 1 1 shows condition numbers for some of the curves shown in figure 8. 
25 The description hereinafter presents simulation results of multiple local control 

systems intended for the active minimization of sound transmitted through a plate. The 
systems are analyzed for harmonic disturbances with respect to stability, convergence, 
reduction of transmitted sound power, the distance between actuators and sensors, and 
sensitivity for reverberating environments. The local control systems are compared with 
30 global control systems. Global control systems are those systems in which each of the 
actuators are controlled in dependence on each of the sensor output signals, whereas 
local control systems are those systems in which one or more of the actuators are con- 
trolled by one or more but not all of the sensor output signals. 
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Figure la shows a baffled plate 1, which supports a plurality of actuators 3(n), 
n = 1, N. In figure la 48 actuators 3(n) are shown. However, if required any other 
number of actuators 3(n) may be applied. 

Supported by suitable supporting means (not shown), a plurality of sensors 2(m), 
5 m = 1, M, is arranged in front of the plate L In figure la, 221 sensors 2(m) are shown. 
Figure la shows a local control system: each actuator 3(n) is associated with 9 sensors 
2(m), adjacent actuators 3(n) sharing three of the sensors 2(m). Of course, any other 
number than 221 sensors 2(m) may be applied and the actuators may be controlled by 
other numbers of sensors. 
10 In figure la, the actuators 3(n) and the sensors 2(m) are regularly arranged in 

columns and rows at equal distances. However, this is not necessary. 

Figure lb shows a cross section through the arrangement according to figure la 
along line IB-IB. The same reference numbers refer to the same elements. 

The acoustic radiation of primary noise source 4 causes a pressure field pmc inci- 
15 dent on plate 1. 

The mutual distance between two adjacent actuators is d x . The mutual distance 
between two adjacent sensors 2(m) is d $ms . The distance between the actuator plane and 
the sensor plane is d. 

Also shown is a reflective wall 8 which might be present in some embodiments, 
20 as will be explained below. 

The actuators 3(n) are shown to be loudspeakers producing secondary noise p s in 
order to reduce the primary noise p p . The total amount of resulting noise is measured by 
the sensors 2(m) which, preferably, are microphones or other pressure-sensitive devices. 

Figure lc shows a schematic electric diagram of the arrangement used in the 
25 invention. The same reference numbers refer to the same components as in figures la 
and lb. 

The sensors 2(m) produce sensor signals p(m) which are transferred to one or 
more controllers 5b(0> * = h X ••*> X e.g., in the way shown in figure lc. 

Figure lc shows four controllers 5b(i), but there may be any other desired num- 
30 ber. They provide one or more output signals Wip which are transmitted to controllers 
5a(i) of a further set of controllers which directly control the actuators 3(n). The outputs 
Wip of the controllers 5b(i) are also input to a supervising controller 6. 

In some embodiments use of one or more detection sensors 7(r), r = 1, R, may 
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be preferred. These detection sensors provide time-advanced infoimation of the primary 
noise p p to a distribution network 10. The distribution network 10 produces detection 
signals v det (i) for the controllers 5a(i). Both the distribution network 10 and the control- 
lers 5a(i) and 5b(i) may be controlled by the supervising controller 6. 
5 Each of the controllers 5a(i) controls one or more of the actuators 3(n) by means 

of control signals u§. 

The supervising controller 6 may be used for monitoring long-term behaviour of 
the system and for modifying control parameters of the distribution network 10 and the 
controllers 5a(i), 5b(i) in order to ensure overall stability of the system. 

10 It is noted that distribution network 10, controllers 5a(i), 5b(i), and supervising 

controller 6 are shown to be separate units, however, in reality they may be implemented 
by a single control unit performing all required functions. Controllers 5a(i), 5b(i) and 6 
are preferably software driven computer units. However, optionally they may be imple- 
mented using digital circuits. Moreover, they need not be physically separated. They 

15 may be implemented as different functional sections of one single processor. On the 
other hand, some of the functionality of their functions may be implemented on remote 
processors if required. For the case of simplicity, in the description and the claims 
reference will only be made to processors 5a(i), 5b(i), and 6. 

Although figure lc shows a situation in which each controller 5a(i) controls one 

20 actuator 3(n), in the theoretical analysis given below, it will be assumed that each con- 
troller 5a(i) controls K actuators 3(n). 
Analysis 

It is assumed that each of the controllers 5a(i), 5b(i) tries to minimize a cost 
function based on sensor signals local to that controller. The scalar cost functions Jj for 
25 the I controllers 5a(i) are written as 

J i =(W i p) H (W i p) + u l H p i u i , i = 0) 

in which p is an M x 1 vector of sensor signals, Wi is a weighting matrix of dimensions 
30 P x M which provides a selection and weighting of P out of a total of M sensor signals 
used as error inputs for controller 5a(i); Ui is a K x 1 -dimensional control signal for node 
i and B, is a K x K dimensional effort weighting matrix. The sensor signals p result from 
the superposition of primary field contributions p p and the contributions p s due to N 
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actuators. The latter contributions are given by Gu, where u is an N x 1 vector denoting 
the control signals that drive the actuators and G is an M x N matrix of transfer functions 
between control signals and sensor signals. Hence, 

P = p p +Gu (2) 

5 

Each controller 5a(i) drives K actuators, so N = IK. 
Introducing the M x N matrix 

G^G^Vfi^Vfi,] (3) 

withFi = Wi H Wi 

10 and Gj denoting the columns G corresponding to controllers 5a(i) having dimensions 
MxK 

and the N x N block-diagonal matrix P defined by 



0...0 

o p 2 ...o 



(4) 



P = 

0 0„.p, 

a linear system of N equations in u can be formulated: 
15 (G H G + P)u = -G H p p (5) 

The present result explicitly includes the weighting factors for the error sensors. 
To arrive at the solution for u an iterative procedure is implemented in the system, such 
as the procedure described by Elliott et al. [5]. For interpretation of system behaviour the 
reader is referred to [1], 

20 Simulations 

In this section simulation results are given for an active control system intended 
to reduce the noise transmitted through plate 1 . The sensors 2(m) are pressure sensors 
placed in the near-field of the plate 1 . In the example, the actuators 3(n) are loudspeakers 
which are assumed to operate as constant volume velocity (monopole-like) sources. The 

25 plate 1 is assumed to be a 1 mm thick aluminium plate of 60 cm x 80 cm, having a 
modulus 7 x 10 10 Pa, Poisson ratio of 0.3, hysteretic damping i\ = 0.02, and a density of 
2.6 x 10 3 kg m" 3 . The plate 1 is assumed to be simply supported and the incident field 
p mc is a plane wave arriving at a direction a of 60 degrees to the plate normal. The basic 
configuration consists of 6 x 8 - 48 actuators and 13 x 17 = 221 sensors, as shown in 
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Fig. la. 

As opposed to active global control systems which minimize a global quadratic 
error criterion, stability is not guaranteed in multiple local systems. Assuming an itera- 
tive procedure to solve Eq. [5], the system is stable if the real parts of the eigenvalues X„, 

5 n = 1, N of the matrix G H G + p are positive [1]. The effort weighting matrix is taken 
to be the diagonal matrix [J = pi. If the system is unstable for p = 0 the value of p will be 
set equal to - min n ReX. n , which makes the system just stable. Increasing the value of IS 
further would enhance the stability margin and improve the speed of convergence of the 
iterative procedure, but also increase the residual radiated power. The convergence of 

10 some iterative procedures is governed by the ratio of the largest singular value ki to the 
smallest singular value k n [5], i.e. the condition number of the Hessian matrix G G + B 
[6]. 

Simulation methods 

The models describing the vibration of the plate 1 can be found in [7]. The pres- 

15 sure p p and p s were computed with a weak form of a Fourier-type extrapolation tech- 
nique in which singularities were evaluated by analytical integration [8]. In principle, the 
Boundary Element method as described in [9] can also be used but the latter method is 
less efficient for geometries of this and larger size. Formulas for zero extrapolation dis- 
tance which were used can be found in [10]. 

20 Simulation results 

The sound power without control and with control for various configurations are 
shown in Fig. 2. It was found that reductions could be obtained for frequencies for which 
both the mutual distance d sen s between the sensors and the mutual distance d x between 
the actuators were smaller than approximately half of a wavelength. Moreover, the dis- 

25 tance d between the sensors 2(m) and the plate 1 turns out to be an important parameter. 
Larger reductions are obtained if the pressure sensors 2(m) are moved away from the 
plate 1. This distance d can not be made arbitrarily large because of stability issues. The 
point of instability is reached at approximately a quarter of a wavelength from the plate 
if the ratio d/d x is larger than a certain minimum value. If this ratio is smaller than this 

30 value, then the system is stable for all frequencies. 

In figure 2 and the other figures, the following notations are used: 

• 48 x 48: 48 sensors and 48 actuators (global control); 

• 48 x 48, 1 x 1 : 48 sensors and 48 actuators, each actuator being depen- 
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dent on one sensor only; 

• 221 x 48: 221 sensors and 48 actuators (global control); 

• 221 x 48, 9 x 1 : 221 sensors and 48 actuators; each actuator being 

dependent on 9 sensors. 
5 A large distance d might be detrimental for primary signals with short correlation 

lengths. For that purpose it may be useful to add one or more detection sensors 7(r) in 
the near-field of the plate. 

The corresponding condition numbers are shown in Fig. 3. If a positive value of 
p was used to make the system stable then the condition number is not shown. 
1 0 Influence of d on the reduction 

From the previous results it was found that the distance d between actuator plane 
and the sensor plane has a considerable influence on the achievable reduction of radiated 
sound power. It was also found that the distance d determines the frequency above which 
the system has to be stabilized by increasing the value of p. A higher value of P leads to 
1 5 smaller reductions. The distance for instability is reached at approximately a quarter of a 
wavelength. 

Clearly, two contradicting requirements for d have to be satisfied for broadband 
reductions. This is illustrated in Fig. 4, which shows sound power radiated from plate 1 
without control and with local control using a 48 x 48, 1 x 1 system, i.e., using a total of 

20 48 sensors and 48 actuators, 1 sensor and 1 actuator for each independent controller, 
with the distance d between the actuator plane and the sensor plane as parameter. If, at 
any frequency, the system is unstable a positive value for p is used which makes the 
system just stable. If the system is stable p = 0 is used. It can be seen that, for small d, 
reductions are increased by increasing d, particularly at low frequencies. However, the 

25 system has to be stabilized above the frequency where d equals a quarter of a wave- 
length. This stabilization leads to smaller reductions at high frequencies. 

Hence, for broadband applications there might be an optimum value for d if the 
objective is to minimize the total acoustic power within a wide frequency range. It is 
assumed that all frequencies are taken into account for which half of the wavelength is 

30 larger than the actuator spacing d x . For the present configuration, this corresponds to all 
frequencies smaller than f < c/2,d x = 1715 Hz. The latter frequency is indicated by a 
dashed line in Fig. 4. This frequency is the maximum frequency for which an active 
control system using a global error criterion leads to significant reductions of radiated 
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sound power. For the present lxl system, the sensor spacing is identical to the actuator 
spacing. The broadband reductions for various values of d normalized to actuator spac- 
ing d x are shown in figure 5. Indeed, it can be seen that there is a maximum in the reduc- 
tion of broadband radiated sound power, both for constant weighting and for 
5 A-weighting. The maximum reduction is obtained for d x /2 < d < d x . 

Additional factors might influence the optimum for d. In the case of stochastic 
disturbances and no reference sensor 7 in a feedforward link, the delay between the ac- 
tuator and the sensor should be small compared to a characteristic correlation length of 
the disturbance signal. In addition, for smaller d, the condition number ki/k n of the sys- 

1 0 tern is lower, and often, therefore, the convergence of adaptive schemes better. These 
two considerations can lead to an optimum for d which is somewhat smaller than given 
by figure 5. Then, for most systems occurring in practice, the optimum for d is in the 
range 0.1d x <d<d x . 

The results for a 221 x 48, 9 x 1 system, having half the distance between the 

1 5 sensors, are shown in Figs. 6 and 7. Figure 6 shows the sound power radiated from a 
plate in such a system, whereas figure 7 shows the broadband reduction, again for all 
frequencies f < c/2.d x . It can be seen that the maximum reduction which can be obtained 
is similar. The optimum value for d, as obtained from figure 7, is also within the range 
d x /2 < d < d x , although the peak in the reduction is wider than in figure 5. In practice 

20 therefore, the value of d for the 9 x 1 system will often be chosen somewhat smaller than 
for the 1 x 1 system. 

The results for a global control system are shown in figures 8 and 9. The differ- 
ences with the preceding local control systems are mainly in the high-frequency range. 
This leads to larger optimum values for d as well as less pronounced maxima. 

25 Comparing figures 5, 7 and 9 with one another shows that in local control sys- 

tems (figures 5 and 7) the reduction of power RP of the total amount of noise relative to 
the primary noise within a predetermined frequency band shows a peak value between 
1/2 < d/d x < 1 whereas in global control systems (figure 9) no substantial peak value is 
present: then, the requirement seems to be 1/2 < d/d x . Moreover, it follows that the 

30 maximum obtainable reduction in all three cases (figures 5, 7, and 9) is similar. Thus, 
when using local control instead of global control a similar value of reduction is pos- 
sible provided the distance between the plane of actuators 3(n) and the plane of sen- 
sors 2(m) is selected carefully, i.e., in the area of the peak value of RP, preferably 
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such that: 



0.9 X RP m ax < RP < RP n 



Performance in reverberating environment 
5 The performance of the local control system was also investigated for the case 

including reflecting parallel plane 8, The distance of this plane 8 to the actuators was 
taken to be 1 m. The reduction which can be obtained with this configuration is shown in 
figure 10 and the corresponding condition numbers in figure 11. It can be seen that for 
reflection coefficients smaller than or equal to 0.9 the control system remains stable and 
10 leads to reasonable reductions. For a reflection coefficient of 0.99 the possible reduction 
above approximately 500 Hz becomes less than for lower reflection coefficients. 
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List of symbols 



G = M x N matrix of transfer functions between control signals u and sensor 
signals p 

5 Gj = MxK matrix of transfer functions between control signals uj and sensor 
signals p. 

Ji = scalar cost function of controller 5(i); 
i=l,2, 

p p = primary field contributions 
10 p = M x 1 vector denoting the sensor signals p(l), p(2), p(m), p(M) 

u = N x 1 vector denoting the control signals u(l) ? u(2), u(n), u(N), that 

drive the actuators 3(n) 
Ui = K x 1 vector denoting the control signals Ui(l), Uj(2), Uj(k), Ui(K) for 

node i 

15 Wj = weighting matrix of dimensions P x M 

Bi = K x K dimensional effort weighting matrix 

k n = smallest singular value of Hessian matrix G H G + p 

k ( = largest singular value of Hessian matrix G H G + p 
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